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1.  INTRODUCTION 

ARCON  Corporation  is  pleased  to  submit  the  following  Final  Report,  which 
describes  the  research  and  modeling  effort  that  has  been  carried  out  under  the  provisions 
of  contract  #  F19628-99-C-0077.  There  are  sections  summarizing  the  work  we  did  on: 
(1)  electron  transport  modeling;  (2)  CEASE[Dichter  et  al.,  1998]  and  Photovoltaic  Array 
and  Space  Power  (PASP)  [Gtissenhoven  et  al.,  1995,  1997]  Dome  dosimeter  modeling; 
(3)  ionizing  particle  energy  deposition  and  coincidence  event  identification  in  the 
Con^Mtct  Environmental  Anomaly  Sensor  (CEASE);  (4)  ionizing  particle  energy 
deposition  in  the  Higb  Energy  Particle  (HEP)  instrument  [Redus,  7PP5];  (5)  proton 
transit  modeling;  and  (6)  simulation  studies  of  grazing  angle  proton  scattering. 

The  overall  objectives  of  the  work  that  we  have  undertaken  are  to:  (1)  perform 
computer  simulations  of  charged  particle  transport,  energy  and  charge  dqmsition  in 
satellite-home  instrumentation  used  in  research  efforts  of  the  Air  Force  Research 
Laboratory,  Space  Weather  Center  of  Excellence  (AFRI7VSBXR)  to  detect  and 
characterize  (by  type,  energy,  intensity,  etc.)  particles  associated  with  ionizing  radiation 
in  space;  (2)  transfer  this  simulation  c^ability  to  AFRiyVSBXR  and  provide  advice  to 
Air  Force  researchers  on  its  use;  (3)  participate  in  the  AFRL  research  effort  by  providing 
the  particle  transport  simulation  component  of  the  research;  and  (4)  enhance  the 
ci^abilities  of  existing  Monte  Carb  transit  programs  by  modifying  these  programs  to 
satisfy  specific  requirements  of  the  AFRL  research  effort.  We  believe  that  these  stated 
goals  are  well  served  by  the  work  reported  bebw. 

Our  rqmrt  is  accon^anied  by  two  pq)ers,  a  rqnint  and  a  pre-prmt,  in  which  the 
author  of  this  report  parHc^ated  as  co-author,  that  were  published  in  IEEE  Transactions 
on  Nuclear  Science  during  the  performance  period  of  fids  contract,  which  was  from  10 
August  1999  until  3 1  August  2003.  The  first  of  these  [Brautigam  et  al.,  2001]  deals  with 
the  solar  cycle  variation  of  outa*  belt  electron  dose  at  low-earth  orbit.  Our  participation  is 
this  research  effort  consisted  of  a  Monte  Carb  analysis  of  the  relative  response  of  the 
CEASE  and  PASP  dosimeters  that  were  used  to  perform  the  electron  dose  measurements. 
This  pq)er  is  included  as  Appendix  A.  The  other  pesper  deals  with  an  investigation  into 
the  causes  of  the  proton  damage  to  one  of  tlK  Chandra  X-Ray  telescope  focal  plane 
detectors.  This  pq)er  demonstrates  how  use  of  standard  proton  transport  Monte  Carb 
calculations  of  grazmg  angle  proton  scattering  produces  significant  under  prediction  of 
the  proton  flux  at  the  tebscope  focal  plane.  This  paper  [Dichter  and  Woolf,  2003] 
conq>rises  .^)pendix  B. 
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2.  ELECTRON  TRANSPORT  MODELING 

The  modeling  of  electron  transport,  energy  and  charge  deposition  was 
accon^lished  using  two  Monte  Carlo  codes,  the  ITS  3.0  Code  System  [Halbleib  et  al., 
1992]  and  the  multi-particle  transport  code  MCNPX  [Waters,  1999].  The  first  of  these, 
ITS  (Integrated  TIGER  Series  of  Co\q)led  Electron/Photon  Monte  Carlo  Codes) 
incorporates  detailed  physical  models  for  electron  and  photon  scattering  and  transport, 
secondary  electron  production,  bremsstrahlung  production,  straggling  and  knock-on 
electron  production.  The  second  code,  MCNPX,  incorporates  all  the  advantages  of  the 
Los  Alamos  MCNP  code  system  [Briesmeister,  1997],  such  as  problem  geometry 
q)ecification  and  code  documentation,  vhile  possessing  the  capability  of  transporting  34 
particle  types  (including  antparticles).  This  code  was  acquired  primarily  for  the  purpose 
of  simulating  proton  tranport,  but  also  proved  to  be  an  effective  tool  for  electrons. 

2.1  Electron  Transport  Simulations  for  the  CEASE  Telescope 

Electron  transport  calculations  were  made  with  the  ITS-ACCEPT  code  for  the 
CEASE  Telescope.  The  ACCEPT  code  geometry  modeling  technique  is  based  on  the 
combinatorial  geometry  method  in  which  the  geometry  model  is  constructed  using  a 
combination  of  body  stupes  offered  by  a  collection  of  basic  body  types:  arbitrarily 
oriented  box(BOX);  rectangular  parallelepped(RPP);  sphere(SPH);  right  circular(RCC) 
and  ellptical  cylinders(REC);  elIpsoid(ELL);  truncated  cone(TO.C);  wedge(^^D); 
arbitrary  polyhedror](ARB).  For  the  CEASE  telescope  without  the  frame  and  case,  a 
total  of 205  todies  and  205  cells  were  defined.  The  ad^ion  of  the  frame  and  case  to  the 
telescope  descrption  increased  the  total  numbers  of  todies  and  cells  to  259.  A  typical 
irput  file  [Woolf,  2000]  sets  ip  a  25000  case  history  ACCEPT  run  for  a  9.9  MeV  electron 
beam  normally  incident  on  the  top  of  the  aluminum  case.  The  beam  cross  sectional  area 
is  a  disk  of  radius  1.4  cm  with  its  center  located  on  the  CEASE  axis  of  symmetry. 

Two  sets  of  ACCEPT  runs  were  made  fra*  the  CEASE  telescope:  telescope 
widiout  and  with  the  frame  and  case.  For  each  configuration,  the  runs  were  made  for  two 
beam  energies,  9.9  and  6.0  MeV,  with  three  incident  obliquities.  The  details  of  these  runs 
and  a  sumnuiryofthe  results  obtained  are  Tqmrted  in  Woolf  [2000].  The  purpose  of  this 
investigation  was  to  obtain  an  estimate  of  the  extent  to  >^ch  the  structure  surrounding 
the  telescope  would  affect  the  resulting  energy  and  charge  deposition  within  the 
telescope. 

Anodier  geometry  model  of  the  CEASE  Telescope  was  also  made  for  use  with 
MCNPX,  and  some  electron  transport  calculations  were  made  paralleling  those  discussed 
above.  TTie  ITS-ACCEPT  and  MCNPX  geometry  models  were  equivalent  with  respect 
to  every  detail. 

2.1.1  CYLTRAN  Simulation  for  the  CEASE  Telescope 

The  CYLTRAN  [Halbleib  et  al.,  1992]  Monte  Carlo  code  is  die  2-dimensional 
cylindrical  geometry  conponent  of  the  ITS  code  series.  It  is  rpplicable  only  in  situations 
v^ere  the  problem  geometry  can  be  ^proximated  assuming  total  cylindrical  symmetry. 
This  symmetry  condition  was  not  strictly  true  for  the  CEASE  telescope,  however 
considerable  economy  of  conputational  effort  was  achieved  because  the  CYLTRAN 
code  which  requires  for  less  conputer  time  than  ACCEPT  if  the  problem  geometry  can 
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be  adequately  {q>proxiinated  in  this  context.  Nearly  equivalent  results,  reported  in  Woolf 
[2000]  were  obt^ed  with  both  codes. 

2.2  Electron  Energy  Deposition  Calculations  in  Silicon  Wafers 

Transport  calculations  were  made  for  4  MeV  and  6  MeV  electrons  incident  on 
rectangular  silicon  dosimeter  wafers  (0.05cm  x  0.9cm  x  0.9cm).  The  Monte  Carlo  runs 
were  made  using  200,000  case  histories  for  all  source  configurations  and  both  source 
energies.  The  sources,  all  located  on  the  wafer  surfeces,  consisted  o£  normally  incident 
electron  beams  (or  point  sources);  point  isotropic  sources;  disk  sources-normal  incidence; 
disk  sources-isotropic  incidence.  The  source  configurations  are  illustrated  in  Table  1. 
Defeult  values  for  the  electron  tow  energy  cut-ofif  were  used  with  both  ACCEPT 
(O.OSEsourec)  and  MCNPX  (1.0  keV).  The  defeult  value  for  the  MCNPX  cut-off  energy 
was  set  much  tower  than  the  ACCEPT  value,  restilting  in  run  times  for  MCNPX  (~1  hr) 
exceeding  those  for  ACCEPT  (~0.5  min).  When  the  same  electron  cut-off  energy  was 
used  in  MCNPX,  the  run  time  was  foimd  to  be  a  fector  of  1.5  greater  than  that  required 
for  ACCEPT. 

Diq>licate  simulation  runs  were  made  with  both  FTS-ACCEPT  and  MCNPX.  Two 
objectives  were  acconq)lished  with  these  calculations:  (1)  a  set  of  input  files  was 
provided  for  both  simulation  programs  that  could  be  modified  and  used  by  AFRL 
personnel;  and  (2)  the  relative  advantages  and  disadvantages  of  the  ACCEPT  and 
MCNPX  codes  for  electron  transport  were  conq>ared.  The  ii:q)ut  files  were  also 
configured  to  produce  electron  pulse-height  spectra  outputs.  We  provided  interpretation 
of  the  pulse-height  spectra  results  for  both  codes.  San:q)le  energy  deposition  results 
obtained  with  ACCEPT  for  the  6  MeV  sources  are  shown  in  Figures  1  and  2. 
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Table  1. 

Eight  Source  Configurations  for  Transport  Simulation  of  4.0  and  6.0  MeV  Electrons  in  Silicon 

[Woolf.  2002\ 


Electron  Source 
Configuration 


(1)  Single  Beam,  Normal 
Incidence  at  (.45, .45,0.) 


Source  Geometry 


(2)  Disk  Source, 

Rad.  =  0.1cm, 

Normal  Incidence  centered 
at  (.45,.45.0.) 


(3)  Disk  Source, 

Radius  =  0.2cm, 

Normal  Incidence  centered 
at  (.45..45,0.) 


(4)  Disk  Source, 

Radius  =  0.449cm, 

Normal  Incidence  centered 
at  (.45, .45.0.) 


(5)  Point  Isotropic,  45^* 
cone  centered  at 
(.45, .45,0.) 


(6)  Disk  Source, 

Rad.  =  0.1cm, 
isotropic  45°  cones 
centered  at  (.45,. 45,0.) 


(7)  Disk  Source, 

Rad.  =  0.2cm, 
isotropic  45°  cones 
centered  at  (.45,.45,0.) 


(8)  Disk  Source, 

Rad.  =  0.449cm. 
Isotropic  45°  cones 
centered  at  (.45,. 45,0.) 


mMmmmi 

•MKtmr 
ti^pi 


Normal  Beam  at  (.45,  .45, 0.)  Normal  Disk,  R^^O.I  at  (.45,  .45, 0.)  Normal  Disk,  RbO.2  at  (.45,  .45, 0.)  Normal  Disk,  R=:0.449  at  (.45,  .45, 0.) 
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Figure  1.  Total  Energy  Deposition  [MeV]  in  Si  Wafer  (0.05x0.9x0.9cm^)  for  8  Source  Geometries;  Electron  Energy  =  6  MeV 
[Woolf,  2002], 


23  ITS-ACCEPT  Electron/Photon  Transport  Simulations  for  the  HEP  Instrument 

An  ACCEPT  geometry  input  file  for  the  HEP  instrument  was  created  from  the 
manufocturing  drawings  supplied  by  Anqjtek,  Inc.  [Redus,  1998].  The  ACCEPT 
geometry  description  contains  the  same  degree  of  detail  as  the  manufocturing  drawings 
permitted:  (1)  taking  advantage  of  the  speed  and  efiBciency  of  the  ACCEPT  code  for 
performing  coupled  electron-photon  transport  calculations  in  con:q)licated  structures;  and 
(2)  use  of  the  ACCEPT  code  to  confirm  the  validity  of  the  MCNPX  geometry  descr^tion 
by  coinparing  electron  transport  results  obtained  with  the  two  programs.  Unlike 
ACCEPT,  MCNPX  could  also  be  used  to  perform  transport  calculations  for  protons, 
neutrons,  mesons,  and  other  particle  species.  It  was  therefore  important  to  have  the 
ACCEPT  version  of  the  HEP  geometry  to  use  as  an  independent  verification  of  the 
MCNPX  geometry  descr^tion.  The  ACCEPT  input  file  listing  can  be  found  in  Woolf, 
2002,  Appendix  1.  Several  ACCEPT  runs  were  made  to  test  the  robustness  of  this 
geometry  file.  The  piupose  was  to  uncover  “holes”  (errors)  in  the  geometry  specification 
due  to  improper  cell  definition.  Test  runs  of  ACCEPT  were  made  with  the  model 
development  at  graduated  stages  of  conplexity  imtil  the  foil  model  illustrated  in  the 
SABRINA[Fa/i  Riper,  2001]  visualuation  shown  in  Figure  3  passed  all  tests  for 
robustness.  _ 


Figure  3.  Sabrina  [Van  Riper,  2001]  Rendering  of  ITS-ACCEPT  Geometry  Model  of 
HEP  Flight  Sensor. 


BEST  AVAILABLE  COPY 
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2.4  rrS-ACCEPT  Program  Enhancements 

Modifications  to  the  ITS-ACCEPT  Monte  Carlo  program  were  made  to  provide 
the  code  with  additional  cq)abilities  that  would  be  useful  in  the  design  and  evaluation  of 
AFRL  space-bome  radiation  sensors.  These  modifications  fell  into  two  categories:  the 
addition  of  new  source  options;  and  the  addition  of  code  that  permits  the  user  to  view  the 
energy  deposition  contributions  of  individual  electron  tracks.  The  primary  motivation  for 
the  source  option  enhancement  Avas  sl^>plied  by  the  feet  that  the  standard  disk  source 
option  in  ACCEPT  does  not  allow  for  the  i^ecification  of  electron  source  beam  slant 
angles  without  slanting  the  souTce  plane.  Source  electrons  emanating  from  a  plane  with 
off-normal  angles  of  incidence  could  not  originate  at  equidistant  points  from  the  target. 
The  demonstration  of  this  is  given  in  a  set  of  electron  transport  runs  that  were  made  for 
the  aluminum- void-silicon  slab  geometry  shown  in  Figure  4.  Runs  Avith  3.5  MeV 
electron  sources  were  made  for  six  source  disk  slant  angles  An  ACCEPT  input  data  file 
corresponding  with  0  =  45®,Avith  the  radius  R  of  the  disk  source  set  at  1.0  cm.,  is  shoAvn 
in  Figure  5. 

2.4.1  Disk  and  Rectangle  Source  Options 

Code  additions  were  Avritten  and  installed  into  the  ACCEPT  program  to  allow  for 
disk  and  rectangular  spatially  uniform  distributed  sources  (electrons  or  photons).  This 
permits  the  user  to  change  the  slant  angle  of  the  source  beam  direction  Avithout  slanting 
the  source  plane  itself.  The  form  of  the  iiqiut  lines  required  to  activate  these  options  are 
given  in  Section  2.3.1  of  Woolf  2002.  The  code  that  Avas  added  to  ACCEPT  also  prints 
diagnostic  messages  in  the  standard  ITS  format  and  aborts  the  run  if  an  error  in  the  input 
Hata  Ariolates  the  conditions  for  geometric  Aralidity. 

Validation  runs  were  made  for  both  the  rectangle  and  disk  source  cases. 

2.4.2  Individual  Electron  Track  Option 

Another  modification  made  to  the  ACCEPT  program  permits  the  user  to  view  the 
energy  deposition  contributions  of  individual  case  histories  (electron  tracks)  in  as  many 
as  10  problem  geometry  cells.  The  form  of  the  input  line  for  inplementation  of  this 
option  and  sanqile  outputs  are  given  in  Section  2.3.2  of  Woolf  2002. 

The  program  listings  containing  these  modifications  to  ACCEPT  that  incorporate 
the  "RECTANGLE -SOURCE" ,  "CIRCLE -SOURCE"  and  "INDIATIDUAL-HISTS" 

options  are  given  in  Appendix  2  of  Wool^  2002. 
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ZONE  #5  -  SPHERE  =  escape  zone 


Figure  4.  Aluminum  /  Void  /  Silicon  ACCEPT  Problem  Geometry  (not  drawn  to  scale)  With 
Slant  Disk  Source  as  Described  in  the  Input  Data  File  Shown  in  Figure  5.  [Woolf,  2002] 
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Figure  5.  ACCEPT  Iiq)ut  Data  File  for  3.5  MeV  45®  Slant  Disk  Source  Incident  on 
Aluminum  /  Void  /  Silicon  Configuration  Shown  in  Figure  4.  [Woolf,  2002] 

2.5  Dome  Dosimeter  Studies 

ITS-ACCEPT  and  MCNPX  were  used  to  perform  Electron  and  proton  transport 
calculations  for  the  CEASE  DDl  and  DD2  dosimeters  [Dichter  et  al.,  1998]  and  PASP 
Dome  D2  and  D3  dosimeters  [Gttssenhoven  et  al.  1995, 1997].  The  calculations  were  made 
to  assist  in  determining  the  effects  of  differences  in  shield  geometry  on  dose  measurements 
in  the  same  radiation  environment.  With  the  differences  in  shielding  geometry  taken  into 
a  common,  normalized  basis  would  be  provided  for  the  use  of  measurement  data 
obtained  with  the  CEASE  and  PASP  dosimeters  to  study  the  solar  cycle  dependence  of  the 
electron  dose  from  the  outer  radiation  belt. 

The  CEASE  DDl  and  DD2  layout  diagram,  a  flat  rectangular  silicon  diode  resting 
on  an  aluminum  oxide  substrate  mounted  on  an  aluminum  base,  [Brautigam,  2001]  for  the 
is  shown  in  Figure  6.  The  corresponding  ACCEPT  geometry  schematic  is  depicted  in 
Figure  7.  The  dosimeter  is  capped  with  an  aluminum  plate.  The  DDl  and  DD2  assemblies 
differ  only  with  respect  to  dimension. 
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Figure  6.  CEASE  DDl  and  DD2  Dosimeter  Assemblies  [Woolf,  2002]. 


Figure  7.  CEASE  DDl  and  DD2  Dosimeter  Geometry  Schematic  for  ACCEPT  Showing 
Plane  Isotropic  Electron  Source;  Geometry  Zones;  nn  =  Material  Cells.  [Woolf, 

2002] 
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2.6  CEASE  Dosimeter  Models 

Models  for  the  original  CEASE  DDl  and  DD2  geometries  were  obtained  directly 
using  the  "CIRCLE-SOURCE"  option  (see  Section  2.4.2).  Spatially  uniform  isotropic 
source  distributions  of  electrons  on  the  aluminum  disc  cover  (see  Figure  6)  were 
assumed.  A  new  source  option,  "DOME  SOURCE,"  was  written  for  ITS- ACCEPT.  For 
this  option,  it  was  assumed  that  an  isotropic,  inward-directed  electron  source  is  uniformly 
distributed  on  the  surface  of  a  thin  void  hemispherical  shell  of  radius  equal  to  that  of  the 
disc  shield  plate  as  shown  in  Figure  8  The  DDl  and  DD2  dosimeters  were  also  studied 
using  MCNPX  for  a  consistency  check.  The  MCNPX  geometry  schematic  is  shown  in 
Figure  9.  The  run  files  for  both  ACCEPT  and  MCNPX  simulations  are  given  in 
Appendices  5  and  6  of  Woolf,  2002. 
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Figure  9.  CEASE  DDl  and  DD2  Dosimeter  Geometry  Schematic  for  MCNPX  Showing 
Hemispherically  Shaped  Isotropic  Electron  Source;  =  Surfaces;  tin  =  Material  Cells, 


\Woolf,  20021 


2.7  PASP  Dosimeter  Models 

The  ITS-ACCEPT  and  MCNPX  codes  were  also  used  to  model  the  PASP  dome 
dosimeters.  The  geometry  schematics  and  source  configurations  for  PASP  Dome  2  for 
both  ACCEPT  and  MCNPX  shown  in  Figures  10  and  11,  respectively.  While  the  use  of 
both  ITS-ACCEPT  and  MCNPX  may  have  been  redundant,  MCNPX  afforded  the 
advantage  of  a  capability  for  modeling  proton  transport  in  addition  to  electron  transport. 
Electron  transport  simulations  were  run  at  AFRL  using  the  CEASE  and  PASP  dosimeter 
models  described  above  for  several  power  law  electron  energy  spectra  characteristic  of 
the  outer  belt  electron  environment.  The  paper  [Brautigam  et  ai,  2001]  in  which  these 
results  are  reported  comprises  Appendix  1  of  this  report. 
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2.8  CRRES  Dosimeter  Modeling 

In  1992,  a  series  of  Monte  Carlo  calculations  were  made  at  Los  Alamos  National 
Laboratory  [Auchampaugh  and  Cayton.,  1992\  in  \riiich  flux  and  dose  response  functions 
were  obtained  for  electrons,  secondary  photons  and  protons  in  the  CRRES  (Combined 
Release  and  Radiation  Effects  Satellite)  dosimeters.  At  the  time  that  these  calculations 
were  made,  the  most  up-to-date  conqjuter  codes  available  were  MCt^lBriesmeister, 
1986]  for  electron  and  photon  transport  and  LAHET[Prflie/  and  Lichtenstein,  1989.]  for 
the  proton  transport.  Because  these  codes  have  since  been  superceded  by  MCNPX,  we 
repeated  some  of  the  Auchampaugh  and  Cayton  electron/photon  transport  calculations 
using  both  the  MCNPX  and  ACCEPT  codes,  and  in  so  doing,  provided  an  “in-house” 
capability  for  Air  Force  persoimel  to  perform  these  and  other  simulations  of  this  kind. 

The  original  MCNP  geometry  files  from  1992  for  the  CRRES_D1,  D2,  D3  and 
D4  dosimeters,  supphed  by  AFRL  [Brautigam,  2002]  were  fed  into  the  visual  editor 
program  VISED[Carter  and  Schwarz,  1997]  for  interpretation,  visualization  and 
translation  to  the  code  ii^ut  fi)rm  of  MCNPX  and  ITS-ACCEPT.  Acceptance  of  the 
original  MCNP  input  file  by  MCNPX  turned  out  to  be  somewhat  problematic  since  there 
had  been  several  input  format  changes  made  since  1992  to  the  MCNP  code  femily,  and 
there  would  have  been  no  guarantee  that  the  original  file  would  have  been  interpreted  by 
MCNPX  as  originally  intended. 

The  earlier  authors,  Auchampaugh  and  Cayton,  used  a  cosine-weighted  current 
source  (isotropic  flux)  and  made  runs  for  17  electron  source  energies.  We  modified  our 
dome  source  option  in  ITS-ACCEPT  to  accommodate  the  cosine-weighted  source  and 
duplicated  the  original  17  runs.  The  ITS-ACCEPT  code  was  chosen  for  this  task 
because,  on  a  run-for-run  basis,  the  run  time  was  shorter  than  that  required  for  MCNPX. 
The  dose  response  functions  calculated  in  this  manner  were  found  to  be  very  close  to  the 
original  results  for  electron  source  energies  below  1  MeV.  For  source  energies  above  1 
MeV,  the  results  were  2-3  orders  of  magnitude  lower  than  the  low  energy  group  and  were 
approximately  one  order  of  magnitude  higher  than  the  original  values.  Graphs  and  a 
table  summarizing  these  findings  are  shown  in  Figure  12  and  Table  2. 

3.  PROTON  TRANSPORT  MODELING 

Throughout  the  performance  of  this  contract,  the  MCNPX  Monte  Carlo  program 
was  used  to  simulate  proton  transport.  These  investigations  were  facilitated  by  the  fact 
that  the  same  geometry  models  that  were  used  for  the  electron  transport  could  be  used  for 
the  proton  work. 

3.1  HEP  Flight  Sensor  Simulations  of  Proton  Transport  with  MCNPX 

Several  Monte  Carlo  runs  were  made,  using  the  MCNPX  model  of  the  in-flight 
version  of  the  HEP  instrument  (see  Section  2.4),  of  proton  transport  and  energy 
deposition.  A  sample  MCNPX  run  file  is  given  in  Appendix  4  of  Woolf,  2000. 
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1e-1 


Figure  12.  Comparison  of  Flux  and  Dose  Response  Functions  Calculated 
With  ITS-ACCEPT  and  as  Reported  by  Auchampaugh  and  Cayton  [Woolf, 
20031 
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3.2  Simulations  of  Proton  Transport  in  the  CEASE  Telescope  with  MCNPX 

The  model  of  the  CEASE  telescope  tiiat  was  written  for  use  with  MCNPX  (see 
section  2.1)  was  also  used  for  studies  of  proton  transport.  Calculations  of  energy 
deposition  in  the  front  (DFT)  and  back  (DBT)  silicon  wafer  detectors  of  the  CEASE 
telescope  were  made  for  eight  proton  source  energies.  The  source  geometry  assumed  was 
a  disk  soirrce  covering  the  telescope  aperture  with  protons  exiting  the  source  plane  along 
a  direction  parallel  to  the  telescope  axis.  The  MCNPX  energy  depositioii  results  (Fip-e 
13)  for  the  proton  sources  nearly  exactly  matched  earlier  energy  deposition  calculations 
supplied  to  us  by  AF91.{Brautigam,  2000\  (Figure  14). 

Proton  Energy  Deposition 
MCNPX  simulations 


1E-2  0.1  1  ^B+^ 


Edep  in  DFT  (MeV) 

Figure  13.  Energy  Depositions  Due  to  Protons  in  DBT  vs.  DFT  for  the  CEASE 
Telescope  as  Con5)uted  With  the  MCNPX  Simulation  Program  [Woolf.  2002]. 
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Figure  14.  Conqiuted  Energy  Depositions  due  to  Protons  in  DBT  vs.  DFT  for  the 
CEASE  Telescope  [Brautigam,  2000]. 
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Slant  source  proton  transport  runs  were  also  made  for  the  CEASE  telescope  for 
two  9  MeV  proton  disk  sources  (radii  0.025  and  0.05  cm)  with  40  degree  slant  incidence. 
The  pulse  height  distribution  results  for  the  DFT  and  DBT  detectors  are  shown  in 
Appendix  3of  Wool^  2003. 

In  addition  to  the  proton  energy  dqwsition  calculations,  a  program  was  written  to 
identify  proton  coincidence  events  in  the  DFT  and  DBT  detectors  by  analyzing  the 
individual  proton  track  histories  in  the  transport  simulation.  The  program,  “count.F,” 
written  for  this  purpose,  performed  an  audit  on  the  MCNPX  track  file  *'ptrac”  nd 
identified  the  history  numbers  in  which  the  same  proton  track  entered  both  DFT  and 
DBT.  This  enabled  die  investigator  to  repeat  die  calculation  to  examine  the  details  of  the 
identified  tracks.  Sample  output  from  “count.F’  and  the  program  listing  may  be  found  in 
Figure  12  and  .^ipen^  3  of  Wool^  2002,  reflectively. 

33  Grazing  Angle  Proton  Scattering  Calculations 

An  investigation  was  launched  to  examine  the  treatment  by  standard  Monte  Carlo 
programs  of  grazing  angle  scatter  of  protons  from  material  surfoces.  This  work  was 
motivated  by  the  discovery  of  extensive  radiation  damage  to  one  of  the  cameras  of  the 
Earth-orbiting  Chandra  X-Ray  Telescope.  The  cause  of  the  damage  to  the  C3iandra 
telescope  was  traced  to  the  effect  of  energetic  protons,  traveling  on  trajectories  nearly 
parallel  to  the  telecope’s  mirrors.  Grazing  angle  scattering  off  the  mirror  surfoces  led  to 
the  transmission  of  the  incident  protons  to  the  focal  plane  and  the  resulting  energy 
deposition  in  the  focal  plane  detectors  by  the  transmitted  protons.  The  extent  of  this 
damage  turned  out  to  be  orders  of  magnitude  greater  than  had  been  eigiected.  In  an 
attenfit  to  pin  down  the  cause  of  the  damage,  intensive  simulation  studies  had  been  made 
with  the  GEANT4  [Giani  et  al.,  2003;  Nartcdlo  et  al.,  2001]  Monte  Carlo  code  using 
models  of  the  Chandra  and  the  European  XMM  X-Ray  telescopes.  In  our  investigation 
we  demonstrated  that  energetic  proton  scattmng  models  used  in  the  conventional,  well- 
established  Monte  Carlo  simulation  codes,  such  as  GEANT4  and  MCNPX,  which  are 
ftpropriate  for  particles  traversing  thin  foils  and  tranfmrt  inside  bulk  material,  are  not 
^id  for  predicting  the  scattaring  of  protons  with  grazing  angles  of  incidence.  To 
accomplish  the  demonstration,  we  used  MCNPX  to  treat  energetic  proton  scattering. 

A  series  of  proton  tranf)ort  calculations  were  made  with  MCNPX,  in  which  the 
physics  of  proton  scattering  is  treated  in  a  similar  way  as  is  done  in  Geant4.  These  runs 
were  made  for  several  proton  source  beam  energies,  ranging  fi»m  0.05  MeV  to  1.0  MeV 
incident  on  iridium  and  aluminum  slabs  at  grazing  angles  of  incidence  (O.F,  0.5‘’,1.0‘’). 
The  emergent,  or  reflected  proton  energy-angle  spectra  (see  Figure  14)  were  examined 
thorough,  and  it  was  found  that  these  results  did  not  exhibit  the  specular  reflection 
[Firsav,  1967]  and  bw  loss  of  energy  that  would  have  been  necessary  to  produce  the 
observed  radiation  damage  to  the  Chandra  telescope  camera.  Moreover,  e^qieriments  had 
been  performed  {Winter  et  al.,  1997;  I^andzelter  et  al.,  1992]  in  which  it  had  been 
reported  that  the  grazing  mcidence  protons  never  actually  entered  the  scattering  material, 
as  assumed  by  standard  Monte  Carlo  codes,  but  in  foct  had  essentially  skimmed  the 
surfoce.  Standard  Monte  Carb  programs,  such  as  MCNPX  and  Geant4  treated  the 
scattering  and  energy  bss  of  grazing  incidence  protons  m  the  same  way  as  those  entering 
at  wide  angles.  The  results  of  our  findings  were  described  in  the  pf>er  entitled,  “Grazmg 
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Angle  Proton  Scattering:  Effects  on  Chandra  and  XMM-Newton  X-Ray  Telescopes”,  by 
Bronislaw  K.  Dichter  and  Stanley  [Dichter  and  Woolf,  2003], 
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Figure  15.  Sample  Plots  of  Probability  Density 
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of  Emergent  Protons  from 


Iridium  with  Exit  Polar  An^e  6^, ,  for  Incident  Angle  =  90.5°,  for  Source  Energy  Eu,  =  250  keV  and 
Exit  Energy  Bins  Eou,  =  0-12.5, 12.5-25.0, 25.0-37.5, 37.5  50.0, 50.0-62.5  keV  [Woolf,  2003], 
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4.  SUMMARY 

During  the  period  covered  by  this  contract,  technical  progress  was  achieved  on 
several  fronts  related  to  the  research,  evaluation  and  adaptation  of  particle  transport 
simulation  programs  for  modeling  the  detection  and  measurement  of  space  radiation  by 
space-borne  sensors.  This  work  included  the  construction  of  realistic  flight  sensor 
conputer  models  for  in^lementation  with  several  Monte  Carlo  simulation  programs,  the 
performance  of  several  particle  transport  and  energy  deposition  calculations,  simulation 
of  space-borne  dosimeters,  and  studies  of  the  scattering  of  grazing  incidence  protons 
from  surfaces  of  material  constituents  of  space-borne  X-ray  telescopes.  Details  of  this 
work  are  provided  in  the  technical  sections  of  this  report  and  in  three  interim  scientific 
rqwjrts  published  in  December  2000,  January  2002,  and  February  2003,  and  two 
technical  psqiers,  one  published  in  December  2001  and  one  to  be  published  in  December 
2003. 

We  have  made  all  of  the  con^uter  codes  of  various  authorshp,  including  our  own, 
available  to  researchers  at  the  Air  Force  Research  Laboratory.  We  have  provided  somce 
code,  where  applicable,  and  various  materials,  including  sample  files  and  output  data, 
relating  to  the  use  of  these  codes,  and  have  also  engaged  in  discussions  on  their  use.  We 
also  n^de  modifications  to  some  of  these  Monte  Carto  transport  codes  to  render  them 
more  expropriate  for  applications  specific  to  the  AFRL  Space  Vehicles  Division. 
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Appendix  A 

The  paper  by  Brautigam  et  al.  [2001]  that  conpises  Appendix  A  contains  Monte  Carlo 
simulation  results  that  were  summarized  in  Section  2.7  of  this  report. 
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APPENDIX  A 


Solar  Cycle  Variation  of  Outer  Belt  Electron 
Dose  at  Low  Earth  Orbit 


D.  H.  Brautigam,  B.  K.  Dichter,  K.  P.  Ray,  W.  R.  Turnbull,  D.  Madden,  A.  Ling,  E.  Holeman, 

R.  H.  Redus,  and  S.  Woolf 


Abstract^  The  solar  cycle  dependence  of  dose  rates  in  the  low 
altitude  *homs*  of  the  outer  zone  electron  belt  is  examined  using 
data  from  TSX5/C£AS£  (solar  maximum  epoch)  and 
AP£X/PASP+  (solar  minimum  epoch).  It  is  found,  after 
normalizing  the  dose  rates  to  account  for  the  different  shielding 
geometries,  that  the  ratio  of  PASP  to  C£ASE  dose  rates  is  -4  for 
L>4.  This  is  contrary  to  the  equal  dose  rates  predicted  by  the 
NASA  A£8M1N  and  A£8MAX  models. 

I.  INTRODUCTION 

EASE,  the  US  Air  Force  Compact  Environmental 
Anomaly  Sensor  [1],  was  launched  into  low  earth  orbit 
on  the  Tri-Service  Experiment-5  (TSX5)  on  7  June  2000  and 
began  returning  data  the  following  day.  The  experiment  was 
integrated  and  flown  by  the  DoD  Space  Test  Program,  and  is 
expected  to  operate  through  October  2001.  Its  13  months  of 
operation  to  date  has  spanned  a  period  of  solar  maximum 
peak  activity.  The  TSX5  low  earth  orbit  is  similar  to  that  of 
the  APEX  spacecraft  on  which  the  PASP+  dosimeter  [2],[3] 
operated  during  the  two  years  (1994  to  1996)  leading  up  to 
solar  minimum.  These  two  data  sets  provide  the  opportunity 
to  study  the  solar  cycle  dependence  of  flie  outer  radiation  belt 
electron  dose. 

It  is  well  known  that  the  highly  variable  outer  radiation  belt 
MeV  electron  population  contributes  to  spacecraft  system 
degradation  through  deep  dielectric  charging  and  radiation 
dose  [4],[5].  This  particle  population  also  represents  a 
radiation  hazard  to  astronauts,  particularly  if  extravdhcular 
activities  are  scheduled  during  magnetically  disturbed 
conditions  when  the  outer  radiation  belt  is  inflated  with 
intense  fluxes  [6].  Thus,  there  is  a  great  incentive  to  improve 
understanding  of  the  dynamics  of  this  particle  population,  and 
develop  more  accurate  models  with  which  to  specify  this 
environment. 
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(jeosynchronous  spacecraft  provide  the  necessary  long¬ 
term  continuous  database  for  studying  solar  cycle  variations 
of  the  outermost  boundary  of  the  outer  belt.  Such  studies  have 
found  that  during  the  two  to  three  years  preceding  solar 
minimum,  the  MeV  electrons  at  geosynchronous  altitude 
exhibit  episodic  enhancements  with  a  27-day  periodicity 
attributed  to  the  well  structured  high  speed  solar  wind  streams 
(HSSWS)  which  flow  from  coronal  holes  near  the  solar 
ecliptic  plane  during  this  phase  of  the  solar  cycle  [7], [8], [9]. 
This  periodic  behavior  has  also  been  observed  in  dose 
measurements  at  low  altitude  polar  orbit  [10],[ll],[12].  This 
pre-solar  minimum  enhanced  activity  results  in  yearly- 
averaged  geosynchronous  fluxes  of  >3  MeV  electrons  which 
can  be  --8  times  higher  than  the  corresponding  averages 
around  solar  maximum  [9]. 

The  trapped  electron  environment  has  been  characterized 
by  long-term  average  models  (NASA  solar  minimum 
AE8MIN  and  maximum  AE8MAX  [13])  and  magnetic 
activity  dependent  models  (APEXRAD  [3]  ,  CRRESRAD 
[14],  and  CIRRESELE  [15]).  The  solar  cycle  phase 
dependence  of  geosynchronous  electrons  referred  to  earlier  is 
inconsistent  with  the  NASA  AES  electron  models  which  show 
that  the  geosynchronous  (I>=6.6)  integral  flux  of  >3  MeV 
electrons  during  solar  maximum  and  solar  minimum  are 
equivalent  [13].  Other  deficiencies  of  currently  available 
models  have  been  previously  discussed  [16],[17].  Spacecraft 
design  engineers  are  interested  in  not  just  long-term  average 
radiation  belt  models  (either  current  or  some  future  improved 
version),  but  models  which  give  the  probability  of 
encountering  various  levels  of  radiation  as  a  function  of  solar 
cycle  phase. 

This  study  examines  the  solar  cycle  dependence  of  dose  in 
the  low  altitude,  high  latitude  boms'  of  the  outer  zone 
electron  belt.  Both  mission-average  dose  and  probability 
distribution  models  are  created  for  the  APEX  and  TSX5 
epochs  which  encompass  solar  minimum  and  maximum, 
respectively.  These  models  show  that  for  the  epochs  studied, 
the  outer  zone  average  dose  rate  from  electrons  is  a  fector  of 
-4  higher  (for  L>4),  and  peak  at  hi^er  L,  during  solar 
minimum  than  during  solar  maximum.  Section  II  describes 
the  relevant  details  of  the  spacecraft  orbits  and  dosimeters, 
and  section  III  discusses  the  construction  of  the  models. 
Section  IV  discusses  the  results,  and  section  V  concludes  with 
a  summary. 
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n.  SATELLITE  ORBITS  AND  DOSIMETERS 

The  APEX  satellite  was  launched  on  5  August  1994  into  a 
70®  inclination  orbit  with  a  perigee  of  362  km,  an  apogee  of 
2544  km,  and  an  orbital  period  of  115  minutes.  The  mission 
had  some  downtime,  but  remained  active  through  May  1996 
(solar  minimum),  with  the  PASP+  dosimeter  returning 
approximately  12  months  of  data.  The  TSX5  satellite  was 
launched  on  7  June  2000  into  a  69®  inclination  orbit  with  a 
perigee  of  410  km,  an  apogee  of  1710  km,  and  an  orbital 
period  of  105  minutes.  The  CEASE  dosimeter  has  been 
returning  data  continuously  (with  occasional  minor  data  gaps) 
to  date,  resulting  in  more  than  12  months  of  data  for  analysis. 
The  solar  cycle  phase  of  the  APEX  and  TSX5  satellite  epochs 
are  shown  in  Fig.  1.  The  APEX  epoch  spans  the  -^2  years 
leading  into  solar  minimum  (May  1996),  and  the  TSX5  epoch 
includes  -^1  year  about  solar  maximum. 

Both  the  PASP+  and  CEASE  dosimeters  are  constructed 
with  an  electrically  active  silicon  diode  sensor  behind  an 
aluminum  (Al)  shield  of  a  prescribed  thickness  which 
determines  the  threshold  energy  of  particles  which  can 
penetrate  the  shield  and  deposit  energy  in  the  silicon  sensor. 
The  PASP+  dosimeter  instrument  includes  four  individual 
dosimeter  heads  (D1  JD2,D3,D4)  providing  measurements 
tirom  a  range  of  particle  energies  [2], [3].  The  D2  dosimeter 
head  includes  two  different  sized  silicon  diode  sensors 
beneath  an  aluminum  dome  shield,  D2B  being  the  larger  of 
the  two.  The  CEASE  instrument  includes  two  dosimeter 
heads  DDl  (thin  shield)  and  DD2  (thick  shield)  [1].  Table  1 
summarizes  the  main  characteristics  of  the  PASP+  D2B 
dosimeter  and  the  CEASE  DDl  dosimeter  (which  for 
simplicity  will  be  leforred  to  from  hereon  as  PASP  and 
CEASE,  respectively)  which  are  relevant  to  this  study. 
Although  the  Al  shield  thickness  fr)r  the  two  dosimeters  is 
comparable,  the  geometry  is  different;  PASP  is  constructed 
with  a  hemispherical  dome  shield  and  CEASE  vyith  a  planar 
shield.  Particles  following  straight  line  paths  towards  the 
sensor  from  any  point  within  the  two  pi  field  of  view  outside 
the  PASP  dome  will  pass  through  equal  pafolengths  of  Al 
shield.  Although  the  planar  dosimeter  does  not  exhibit  this 
spherical  symmetry  it  represents  a  more  realistic  geometry  for 
simulating  spacecraft  shielding  designs.  The  PASP  and 
CEASE  silicon  diodes  are  of  comparable  thickness,  yet  the 
CEASE  diode  has  a  much  greater  planar  dimension.  The  Low 
Linear  Energy  Transfer  (LoLET)  channel  is  defined  by  its 
energy  deposition  threshold  settings  (0.05  to  1.0  MeV  for 
PASP  and  0.05  to  0.85  MeV  for  CEASE).  For  ^80  mils  of  Al 
shielding,  electrons  above  -1  MeV  and  protons  above  78 
MeV  (for  CEASE)  to  1 15  MeV  (for  PASP)  can  penetrate  the 
shielding  and  deposit  the  required  energy  for  the  LoLET 
channel.  Although  the  >1  MeV  electrons  are  the  dominant 
source  of  LoLET  dose  in  the  outer  radiation  belt,  the  sporadic 
solar  energetic  proton  events  (most  fi*equent  during  solar 
maximum)  may  contribute  a  relatively  significant  short  term 
LoLET  dose  over  the  course  of  a  2  to  3  day  event.  These 


events  have  been  identified  and  are  excluded  firom  the 
database  for  this  study. 

Dose  measurements  are  dependent  on  shielding 
characteristics  (material,  thickness,  and  geometry)  and  thus  do 
not  provide  a  direct  measure  of  the  environment  as  does  a 
fiux  measurement  which  discriminates  particle  species  and 
energy.  In  order  to  substantiate  any  claim  about  the  variable 
electron  environment  based  on  the  comparison  between  two 
different  dosimeter  measurements  on  separate  spacecraft,  it  is 
important  to  'normalize'  their  dose  values  to  some  reference. 
The  relative  response  of  the  two  dosimeters  is  investigated 
using  the  Integrated  Tiger  Series  (ITS)  electron  transport 
code  [18].  In  this  simulation,  an  isotropic  power  law  spectrum 
(I  to  10  MeV)  is  uniformly  emitted  firom  a  virtual 
hemispherical  surfece  centered  on  the  dosimeter  sensor. 
Simulations  are  performed  for  a  range  of  power  law  indices 
(N=2  to  6)  vMch  are  representative  of  the  outer  belt  electron 
environment.  The  results  show  tiiat  for  the  hardest  incident 
flux  spectrum  (N=2)  the  dose  ratio  for  the  hemispheric  to 
planar  dosimeter  is  2.0,  and  increases  monotonically  to  2.4 
for  the  softest  spectrum  (N=6),  for  an  average  ratio  of  2.2. 
Because  the  planar  dosimeter  better  represents  spacecraft 
surfece  geometry,  and  because  a  number  of  CEASE  units  are 
planned  for  future  flights,  CEASE  is  taken  as  the  reference 
standard  here,  and  the  PASP  (hemispheric)  dose  is 
normalized  to  CEASE  dose  by  dividing  by  a  fector  of  2.2. 

Another  issue  to  consider  when  comparing  particle  data 
fixnn  different  instruments,  even  fix)m  the  same  spacecraft,  is 
their  field  of  view  (FOV)  relative  to  the  geomagnetic  field. 
This  determines  the  particle  pitch  angle  distribution  being 
sampled  (the  pitch  angle  is  the  angle  between  the  particle 
velocity  and  the  local  geomagnetic  field  vectors).  A 
collimated  particle  telescope  pointing  in  a  direction 
orthogonal  to  the  geomagnetic  field  can  see  orders  of 
magnitude  higher  particle  fluxes  (at  90®  pitch  angle)  than  an 
identical  instrument  pointing  parallel  (anti-parallel)  to  the 
field  line  at  0®  (180®)  pitch  angle.  The  APEX  and  TSX5 
spacecraft  are  both  3 -axes  stabilized.  The  APEX  attitude 
control  system  maintained  a  solar  inertial  reference  with 
PASP  pointing  nearly  parallel  (within  5®)  to  the  ecliptic  plane 
and  orthogonal  to  the  Earth-Sun  line.  The  TSX5  attitude 
control  system  maintains  a  geocentric  reference  with  CEASE 
pointing  in  a  direction  orthogonal  to  the  spacecraft  velocity 
vector  and  pointing  30®  from  the  zenith  direction.  Assuming  a 
collimated  FOV  direction  normal  to  the  dosimeter  face  (in 
actuality,  each  dosimeter  has  a  full  two  pi  steradian  FOV),  a 
pitch  angle  histogram  was  created  and  a  mission  average 
computed  for  each  dosimeter.  With  the  0-180®  pitch  angle 
distribution  collapsed  to  a  0-90®  distribution,  the  average 
values  were  found  to  be  35®  (CEASE)  and  50®  (PASP)  over 
the  range  ftom  L  equal  3  to  6.5.  This  marginal  difference  in 
average  pitch  angle  represents  a  negligible  effect,  particularly 
given  the  feet  that  the  two  dosimeters  are  not  collimated  but 
have  a  full  two  pi  steradian  FOV.  The  difference  in  PASP  and 
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CEASE  orientation  is  thus  assumed  to  be  an  insignificant 
factor  in  comparing  tiie  daily  average  dose  rates. 

m.  ANALYSIS 

Since  radiation  belt  particle  motion  is  tied  to  the 
geomagnetic  field,  magnetic  coordinates  (L3/B0)  are  the 
logical  choice  for  ordering  the  dose  data.  The  coordinate  ‘L* 
labels  a  magnetic  field  line  by  the  radial  distance  (in  Earth 
radii)  at  which  it  crosses  the  magnetic  equator.  The 
coordinate  'B/Bo’  is  the  ratio  between  the  magnetic  field  B  at 
some  arbitrary  point  along  a  field  line  and  the  equatorial 
(minimum)  value  Bo  along  that  same  field  line.  The  CEASE 
and  PASP  dose  databases  are  binned  and  averaged  in 
(L,B/Bo)  in  a  maimer  identical  to  the  published  APEXRAD 
model  (see  [2],[3]  for  details).  The  resulting  mission  averaged 
APEX/PASP  and  TSX5/CEASE  dose  rate  maps,  when 
displayed  on  the  northern  quadrant  of  a  dipolar  grid  (Fig.  2), 
aid  in  visualizing  the  relevant  spatial  regions.  The  inner  zone 
is  seen  as  the  darker  color  contours  extending  firom  the 
equator  inside  I>=1.5,  while  the  outer  zone  lioms'  are  seen  as 
darker  color  contours  at  high  latitudes  in  the  region  between 
L=3  to  5.  The  heart  of  the  outer  radiation  belt  is  the  equatorial 
projection  of  the  L=3  to  5  field  lines  and  exhibit  much  higher 
dose  rates  than  those  at  the  low  altitudes  shown  in  Fig.  2. 

Since  the  APEX  and  TSX5  orbits  are  not  identical,  regions 
of  (L,B/Bo)  space  are  sampled  differently  by  the  two 
spacecraft.  Singly  averaging  over  the  entire  range  of  B/Bo  for 
a  fixed  L  would  yield  skewed  results  because  of  the  different 
temporal  coverage  in  (L,B/Bo)  and  the  dose  dependency  on 
B/Bp.  To  best  compare  dose  variations  seen  by  PASP  and 
CEASE,  a  subset  of  foe  complete  (L,B/Bo)  grid  is  determined 
which  provides  anomie  observations  for  the  study,  yet  also 
minimizes  the  spread  in  B/Bp  for  a  fixed  L  bin  width  (0. 1  R^). 
The  bottom  panel  of  Fig.  3  illustrates  foe  increasing  B/Bp  bin 
width  as  a  function  of  L.  The  top  panel  illustrates  foe  range  of 
altitudes  sampled  by  foe  respective  satellites  corresponding  to 
the  specified  L  and  B/Bp  range.  For  L  <  5.5  foe  coverage  in 
altitude  is  quite  constant,  where  the  average  range  for  APEX 
altitudes  is  1355  +A  520  km,  and  for  TSX5, 1270  +/-  350  km. 
Fig.  2  illustrates  roughly  foe  region  of  interest  as  being 
enclosed  between  foe  L  =  3  to  5  field  lines  and  between  the 
1000  to  2000  km  altitude  arcs.  In  what  follows,  it  should  be 
understood  that  within  foe  explicit  L  dependence  of  dose 
distributions  is  an  in^licit  dependenoe  in  B/Bp. 

Subsequent  to  binning  foe  full  resolution  (5s)  dose  data  into 
the  (L,B/Bo)  grid  prescribed  above,  foe  daily  averaged  dose 
rates  for  both  data  sets  are  computed  as  a  function  of  L.  These 
daily  dose  profiles  are  plotted  as  epoch  surveys  in  Fig.  4a 
(APEX/PASP)  and  Fig.  4b  (TSX5/CEASE),  along  vrith  foe 
respective  magnetic  activity  index  Dst  history  for  each  epoch. 
The  PASP  dose  survey,  particularly  up  through  1995.5, 
highlights  foe  periodic  enhancements  of  MeV  electrons 
characteristically  observed  during  the  years  approaching  solar 
minimum.  Dst  is  a  measure  of  magnetic  storm  activity,  with  a 
more  negative  value  indicative  of  greater  activity.  Periods 


exhibiting  the  largest  Dst  excursions  (e.g.,  during  foe  TSX5 
epoch  around  years  2000.55  and  2001.25)  also  show  foe 
largest  increases  in  dose  rates  in  the  slot  region  (dovm  to 
L-2.5). 

The  simplest  way  to  compare  the  dose  rates  from  foe  two 
different  epochs  is  to  construct  a  single  average  dose  model 
as  an  explicit  function  of  L  from  each  mission.  These  average 
dose  profiles  are  presented  in  Fig.  5.  The  average  PASP 
(solar  minimum  epoch)  dose  rate  is  higher  and  peaks  at  a 
higher  L  value  than  that  for  CEASE  (solar  maximum  epoch). 
The  PASP  peak  dose  rate  is  l.lxlO*^  rd/s  (at  L=4.1),  while 
foe  CEASE  peak  dose  rate  is  3.2x10’^  rd/s  (at  L=3.8).  The 
ratio  of  PASP  to  CEASE  dose  rate  decreases  from  --5  for 
L-4.5  to  --3  at  L-6.5.  The  PASP  dose  rate  peak  is  also 
somewhat  narrower  than  the  CEASE  peak.  For  comparison, 
average  dose  rates  at  several  (L,B/Bo)  points  are  confuted 
using  spectra  from  the  NASA  AES  models  (the  AES  models 
were  downloaded  fi'om  the  web  ]3ages  of  NASA’s  National 
Space  Science  Data  Center).  For  a  given  L,  foe  values  of  B/Bo 
used  to  specify  foe  AE8MIN  (AESMAX)  spectra  are  the 
average  B/B^  fi’om  foe  corresponding  APEX  (TSX5) 
ephemeris.  The  ITS  code  is  used  to  simulate  the  transport  of 
the  AE8MIN  and  AESMAX  omnidirectional  flux  spectra  (1 
to  7  MeV)  through  the  PASP  and  CEASE  dosimeter  models, 
respectively.  These  AES  model  predicted  dose  rates  are 
included  in  Fig.  5.  For  L>4,  foe  AESMAX  and  AE8MIN 
models  give  identical  fluxes  for  a  given  coordinate  pair 
(L,B/Bo).  However,  due  to  foe  slight  differences  in  the 
average  B/Bo  for  foe  two  spacecraft  orbits,  foe  AE8MIN  and 
AESMAX  dose  predictions  for  L>4  differ  by  --30%,  foe 
difference  attributed  to  the  variation  of  spectra  with  B/Bo-  For 
L<4.5  the  AESMIN  model  is  less  than  a  fiactor  of  2  higher 
than  foe  observed  PASP  dose.  For  L>4  foe  ratio  of  AESMIN 
to  PASP  increases  steadily  to  --5  near  geosynchronous.  The 
agreement  between  AESMAX  and  CEASE  is  much  less,  with 
foe  ratio  of  AESMAX  to  CEASE  increasingly  steadily  firom 
-6  at  L=3.5  to  ~10  near  geosynchronous. 

In  addition  to  long-term  average  dose  models  which 
obscure  foe  wide  dynamic  variations  of  the  outer  belt 
electrons,  space  system  design  engineers  are  interested  in  foe 
probability  of  exceeding  a  given  dose  rate  at  a  fixed  point  in 
space.  To  construct  these  distributions,  foe  observed  dynamic 
range  of  dose  rates  (spanning  five  orders  of  magnitude,  lO’  to 
I0‘^  rd/s)  is  divided  into  50  equal  logarithmically  spaced  bins. 
For  each  L  bin,  foe  fraction  of  5s  dose  measurements  foiling 
within  each  dose  rate  bin  is  determined  such  that  the  integral 
over  all  dose  rate  bins  sums  to  1.  The  probability  for 
exceeding  a  given  dose  rate  D  is  then  determined  by 
integrating  over  these  dose  rate  bins  from  D  to  Dmax(=10'^ 
rd/s).  These  probability  distributions,  for  both  APEX  and 
TSX5  missions,  for  all  L,  are  shown  in  Fig.  6(a;b).  The  shape 
of  foe  dose  profiles  (in  L)  seen  in  Fig.  5  is  reflected  in  Fig.6. 
To  focilitate  a  direct  comparison  of  foe  probability 
distributions  for  foe  two  epochs,  the  distributions  at  a  single  L 
near  the  dose  peak  (1^4)  are  plotted  together  in  Fig.  7.  From 
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this  plot,  it  can  be  seen  that  for  these  missions,  the  probability 
of  exceeding  the  dose  rate  of  10"^  rd/s  is  --0.38  during  the 
solar  minimum  epoch  (APEX)  and  ~0.10  during  the  solar 
maximum  epoch  (TSX5).  It  should  be  emphasized  that  these 
probabilities  have  been  derived  from  a  specific  range  of 
(L,B/Bo)  space.  For  the  case  here  of  L=4,  B/Bo  ranges 
between  48  and  80,  and  is  not  representative  of  the  dose 
environment  for  L=4  on  the  equator  (B/Bo=l)  where  dose 
rates  are  considerably  larger. 

IV.  DISCUSSION 

One  clear  result  here  is  that  for  I>4,  the  mission  averaged 
PASP  (solar  minimum  epoch)  dose  rates  are  a  foctor  of  -A 
higher  than  those  for  CEASE  (solar  maximum  epoch),  a 
significant  departure  firom  the  nearly  equal  rates  predicted  by 
the  NASA  AES  models.  This  result  is  consistent  with  recent 
studies  which  attribute  periodic  enhancements  of  MeV 
electron  flux  at  geosynchronous  to  the  HSSWS  prominent 
during  the  approach  towards  solar  minimum  [7],[8],[9].  This 
27~day  periodicity  in  enhancements  of  MeV  electrons  is 
evident  in  the  PASP  dose  seen  in  Fig.  4a,  where  the  enhanced 
dose  rates  are  observed  throughout  a  wide  range  of  L 
spanning  the  horns’  of  the  outer  radiation  belt. 

A  second  result  is  that  the  average  dose  profile  for  CEASE 
(solar  maximum  epoch)  peaks  at  a  lower  L  (L=3.8)  and  has  a 
broader  peak  than  does  the  dose  profile  for  PASP  (solar 
maximum  epoch)  which  peaks  at  L=4.l.  This  result  is 
attributed  to  the  greater  number  of  magnetic  storms  occurring 
during  solar  maximum.  The  difference  in  geomagnetic 
activity  for  the  two  epochs  is  apparent  in  the  plots  of  Dst  in 
Fig.  4,  where  maximum  negative  excursion  in  Dst 
corresponds  to  the  severest  magnetic  storm  activity.  Although 
there  were  comparable  numbers  (-^)  of  moderate  storms  (Dst 
=  -100  to  -150  nT)  during  both  epochs,  there  were  8  major 
storms  (Dst  <  -150  nT)  during  the  TSX5  epoch,  with  none 
during  the  APEX  epoch.  Associated  with  these  largest  storms 
are  the  iigections  of  dose  producing  electrons  which  penetrate 
deep  into  the  slot  region  (L~2.5).  The  greater  frequency  of 
these  storms  during  solar  maximum  moves  the  inner  edge  of 
the  average  dose  profile  to  lower  L  and  broadens  the  profile 
peak. 

Future  work  will  investigate  the  possibility  of  extending  the 
spatial  coverage  of  in  situ  dose  measurements  by  developing 
algorithms  for  mapping  (fose  along  field  lines  to  other  renK)te 
regions  and  by  interpolating  between  other  existing 
spacecrafr.  Althougjhi  this  study  dealt  with  the  LoLET  dose 
behind  80  mils  of  A1  corresponding  to  >1  MeV  electrons, 
future  work  will  analyze  a  second  CEASE  channel  for  dose 
behind  250  mils  of  A1  corresponding  to  >3  MeV  electrons. 
The  relative  dose  from  these  two  channels  may  then  be  used 
to  examine  the  spatial-temporal  variability  of  spectral 
hardness  of  this  very  dynamic  trapped  particle  population. 


V.  SUMMARY 

The  significance  of  the  outer  zone  electrons,  both  in  terms 
of  their  effect  on  space  technologies  [4],[5]  as  well  as  human 
safety  [4],[6],  is  well  known.  However,  much  work  is 
required  to  better  characterize  this  population's  variability 
over  solar  cycle  time  scales.  Whereas  long-term  averages 
were  once  the  model  standard,  it  is  becoming  increasingly 
evident  that  more  detailed  information  on  the  extremes  and 
the  variability  of  the  environment  is  advantageous  for  both 
satellite  design  engineers  and  those  responsible  for  planning 
human  activities  in  space. 

This  study  has  produced  results  that  are  consistent  with 
studies  utilizing  geosynchronous  flux  data  which  found  that 
electron  enhancements  minimize  near  solar  maximum  and 
maximize  during  the  approach  to  solar  minimum.  The 
discrepancy  with  the  NASA  AE8  models  underscores  the 
need  to  promote  further  improvement  of  existing  radiation 
models  by  continually  monitoring  the  radiation  belts 
throughout  the  solar  cycle.  This  study  has  also  produced  some 
probabilistic  models  as  a  step  towards  an  alternative  to  foe 
standard  mission  averaged  model.  Although  foe  study's  spatial 
coverage  is  limited  to  low  Earth  orbit,  and  its  temporal 
coverage  to  <3  years  of  an  1 1-year  solar  cycle,  it  does  provide 
a  framework  for  incorporating  foture  data  sets.  More  CEASE 
units  are  scheduled  to  fly  in  foe  future,  and  these  data  will 
continue  to  expand  foe  current  spatial-temporal  coverage 
leading  to  improved  modeling  capabilities. 
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TABLE  1 

DOSIMETER  CHARACTERISTICS 


Dosimeter 

Shield  (Al) 

Sensor  (Si) 

LoLET  (MeV) 

Sensor 

head 

Thick  (mil) 
geometry 

Area 

(cm^) 

Depth 

(cm) 

Energy 

deposit 

Electron 

threshold 

PASP 

D2B 

82.5 

dome 

0.051 

0.038 

0.05-1.0 

^  1.0 

CEASE 

DDl 

80.0 

planar 

0.810 

0.050 

0.05-0.85 

^  1.0 

30 


Year 

Fig.  1.  The  solar  cycle  is  depicted  here  in  terms  of  die  smoothed  sunspot 
number.  ITie  dash^  curve  beginning  in  niid-2001  is  the  predicted  value. 
The  mission  time  intervals  for  APEX  (-2  years  leading  into  solar 
minimum)  and  TSX5  ('-I  year  at  solar  maximum)  are  indicated  by  the 
heavy  solid  lines. 
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Fig.  2.  Mission  average  dose  rates  (rd/s)  in  silicon  (color  coded)  plotted 
in  dipole  (L,B/Bo)  coordinates  for  PASP  (top)  and  CEASE  (bottom).  The 
vertical  and  horizontal  scales  are  in  units  of  Earth  radii.  Dipole  field  lines 
for  L  =  1.5,  2,  3,  4,  and  5  have  been  drawn,  with  the  L=1.5  field  line 
shown  intersecting  the  magnetic  equator.  Dotted  circular  arcs  are  drawn 
at  1(X)0  and  2(X)0  km  altitud 
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Fig.  3.  Dose  data  are  sorted  by  (L,B/Bo)  coordinates.  The  L  bin  width  is 
fixed  at  0.1  Rb,  whereas  the  B/Bo  bin  width  increases  with  L.  In  the 
bottom  panel  the  average  B/Bo(L)  for  the  two  spacecraft  is  plotted 
(heavy  solid  line)  versus  with  the  upper  and  lower  limits  of  the 
B/Bo(L)  bin  width  defined  by  +/-  0.25  of  the  average.  In  the  top  panel 
the  altitude  coverage  (+/-  one  sigma  of  the  average  value)  of  satellites  is 
plotted  versus  L 
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Fig.  4a.  Plotted  is  the  APEX  mission  survey.  The  top  panel  plots  PASP  dose  rate  (rd/s)  in  silicon  (color  coded)  as  a  function  of  L  (implicit  B/Bo  dependence) 
and  year.  The  bottom  panel  plots  the  magnetic  activity  index  Dst  (nT)  as  a  function  of  year. 
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Fig.  4b.  Plotted  is  the  TSX5  mission  survey.  The  top  panel  plots  CEASE  dose  rate  (rd/s)  in  silicon  (color  coded)  as  a  function  of  L  (implicit  B/Bodepradence) 
and  year.  The  bottom  panel  plots  the  magnetic  activity  index  Dst  (nT)  as  a  function  of  year.  The  length  of  the  time  scale  is  the  same  as  that  for  Fig.  4b  to 
focilitate  comparison  of  temporal  variations. 
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Fig.5.  APEX/PASP  (heavy  line)  and  TSX5/CEASE  (thin  line)  mission 
average  dose  rates  (rd/s)  in  silicon  as  a  function  of  L  (implicit  B/Bo 
dependence).  Predicted  dose  rates  ftom  AESMIM  (squares)  and 
AE8MAX  (triangles)  are  included  for  comparison. 


L-Shell 

Fig.  6a.  PASP  (solar  minimum  qxx;h)  probabilistic  dose  model. 
Probability  (color  coded)  for  exceeding  a  some  threshold  dose  rate  (rd/s 
in  silicon)  plotted  as  a  function  of  dose  threshold  and  L. 


L-Shell 

Fig.  6b.  CEASE  (solar  minimum  epoch)  probabilistic  dose  model. 
Probability  (color  coded)  for  exceeding  a  some  threshold  dose  rate  (rd/s 
in  silicon)  plotted  as  a  function  of  dose  threshold  and  L. 


TSX5/CEASE  (+)  APEX/PASP  (❖) 


Fig.  7.  This  plot  is  a  2-D  representation  of  the  information  in  Fig.  6a  and 
6b  for  1^4.0.  It  gives  the  probability  that  the  PASP  (triangle)  and 
CEASE  (cross)  dose  rate  at  L=4  exceeds  some  specified  threshold  level 
(independent  variable). 


BEST  AVAILABLE  COPY 


33 


Appendix  B 

The  paper  by  Dichter  and  Woolf  [2003]  that  coii^)rises  >^)pendix  B  presents  results  of 
the  Monte  Carlo  electron  transport  simulations  that  are  summarized  in  Section  3.3  of  this 
report. 
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APPENDIX  B 


Grazing  Angle  Proton  Scattering:  Effects  on 
Chandra  and  XMM-Newton  X-Ray  Telescopes 


Bronislaw  K.  Dichter  and  Stanley  Woolf  Member,  IEEE 


Abstract — proton  scattering  process  resulted  in  damage  to 
one  of  the  Chandra  X-Ray  teiescope’s  focal  plane  detectors.  In 
this  process,  incident  protons  were  transmitted,  by  scattering 
off  the  telescope  mirrors,  to  the  focal  plane.  We  identify  the 
proton  population  responsible  for  the  damage  and,  using  a 
proper  grazing  angle  formalism,  we  show  that  the  standard 
calculations  of  grazing  angle  scattering  will  significantly  under 
predict  the  expected  proton  flux  at  the  focal  plane. 

I.  Introduction 

EVERAL  weeks  into  the  Chandra  X-Ray  telescope 
mission,  unexpected  damage  was  (4>served  to  one  of  its 
cameras  [1].  The  problem  was  identified  as  radiation 
damage  to  the  front  illuminated  charge  coupled  detectors 
(CCD)  the  Advanced  CCD  Imaging  Spectrometer  (ACIS). 
The  amount  of  damage  was  orders  of  magnitude  larger  than 
was  to  be  expected  this  early  in  the  mission.  This  event 
resulted  in  intensive  study  of  the  transmission  of  protons 
through  the  Chandra  and  the  European  XMM-Newton  x-rsfy 
telescopes  [2]  using  the  Monte  Carlo  computer  code  Geant4 
[3].  Both  of  these  instruments  utilize  grazing  incidence 
mirrors  to  focus  the  x-r^  onto  the  CCD  cameras  in  the 
focal  plane. 

In  this  paper  we  will  identify  the  particle  population 
responsible  for  the  damage  and  show  that  the  transmission 
calculations  of  [2]  use  a  model  of  proton  scattering  beyond 
its  range  of  validity,  leading  to  an  underprediction  of  the 
transmission  probability.  We  will  present  grazing  angle 
spattering  data  and  calculations  and  compare  them  to  results 
computed  using  the  same  methods  as  are  used  in  Geant4. 
Finally,  we  will  discuss  the  effect  that  correct  calculation  of 
grazing  angle  scattering  has  on  the  calculated  fluxes  that 
reach  the  Chandra  and  XMM  focal  planes  and  suggest  a  way 
of  obtaining  more  accurate  results. 

II.  Radiation  Damage  TO  CCD’s  ON  Chandra 

Energy  loss  in  Si  solid  state  detectors,  such  as  CCD’s  &lls 
into  two  classes  ionizing  and  non-ionizing  energy  loss 
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(NIEL).  Ionizing  energy  loss  is  due  to  the  distant  collisions 
of  the  incident  particles  with  the  conduction  band  electrons 
in  the  target  material.  This  lype  of  process  leads  to 
temporary  damage  that  can  be  reversed  with  time  or  by 
annealing  the  material.  The  NIEL  process  consists  of 
collisions  with  the  nuclei  in  the  material  lattice.  Any  such 
collision  that  results  in  the  energy  transfer  greater  than  about 
30  eV  will  knock  the  target  atom  fi-om  its  location  in  the 
lattice.  The  result  of  this  collision,  one  atom  in  an 
interstitial  location  and  one  vacancy  in  the  lattice,  is  called  a 
Frankel  defect  and  is  the  most  common  fype  of  bulk  damage. 
Frankel  defects  act  as  charge  carrier  traps,  removing 
electrons  fi'om  the  charge  collection  process.  Protons  with 
ener^  above  a  few  keV  have  sufficient  energy  to  cause  this 
defect.  The  probabilify  of  causing  the  defect  increases  with 
energy  up  to  about  100  keV  and  then  decreases  with 
increasing  energy  with  an  energy  dependence  of  1/E  [4], 

Gate  Region 


CCD’s  in  the  Chandra  instrument  focal  plane.  Region  A 
corresponds  to  the  first  45  pm  of  Si  fiicing  the  incident 
particles.  Back  illuminated  device  gate  region  is  shielded 
from  transmitted  protons  by  the  upstream  region  A.  Front 
illuminated  devices  do  not  have  such  shielding. 

The  ACIS  configuration  of  its  10  CCD’s  is  shown  in 
Fig.l.  Eight  of  the  CCD’s  are  front  illuminated  (FI),  with 
the  charge  transfer  gate  region  directly  exposed  to  incident 
particles  transmitted  through  the  grazing  optics  of  the 
telescope.  Two  are  back  illuminated  (BI)  and  have  the  body 
of  the  device  shielding  the  gate  region.  ACIS  operates  with 
a  thin  Polyamide  film,  with  a  light  blocking  aluminum 
coating  upstream  of  the  CCD’s.  The  total  mass  density  of 
the  film  and  aluminum  is  sufiScient  to  stop  protons  with  E  < 
80  keV.  During  the  first  few  weeks  of  the  mission,  the  FI 
CCD’s  suffered  a  degradation  of  performance.  The  measure 
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of  the  small  inefiBciency  in  transferring  electrons  from  one 
pixel  to  another  during  the  readout  cycles,  or  charge  transfer 
inefiBdency  (CTI),  was  increasing  fer  more  rapidly  than 
expected  for  ie  FI  CCD’s.  The  CTI  for  the 


« 

K 


Fig  2.  Integral  proton  fluence  and  range  in  Si  plotted  as  a 
function  of  energy.  Protons  models  are  described  in  [S]. 


BI  CCD’s  remained  unchanged.  The  ACIS  CCD  integral 
proton  fluence  computed  for  the  Chandra  orbit  is  shown  in 
Fig.  2.  Our  calculations  assume  that  the  incident  proton 
population  has  an  access  path  to  die  CCD’s  that  does  not 
degrade  its  spectral  shape.  We  will  justify  this  in  Section  4. 
Electron  fluences  are  not  a  concern  because  the  low  energy 
electrons  (E  <  100  keV)are  swept  away  by  the  ACIS  broom 
magnets  and  the  hi^er  energy  electrons  are  both  &r  less 
numerous  and  are  hi^y  inefficient  producers  of  Frankel 
defects.  The  FI  CCD’s  can  be  divided  into  Region  A,  the 
same  length  as  BI  CCD’s,  and  remainder.  Region  B.  Both 
BI  and  FI  regim  A  are  sensitive  to  flux  Fa  =  1.65  x  10^ 
protons/cm^-s  coresponding  to  E>80  keV.  The  region  B, 
however  is  only  sensitive  to  fy  =  2.37  x  10^  protons/cm^-s 
corresponding  to  E>1,900  keV).  Since  the  ratio  fa/fA^  1*4  x 
10*^,  the  damage  to  the  bulk  silicon  in  region  B  is  small 
compared  to  region  A  and  can  be  n^ected.  Taking  into 
account  the  feet  the  BI  CCD’s  did  not  suffer  any  degradation 
it  is  likely  that  the  damage  occurred  near  the  front  of  the 
devices.  This  is  made  more  evident  by  computing  <R>,  the 
average  range  of  protons  in  Si,  wei^ted  the  incident 
particle  flux  fl[E), 


<R>  = 


jR(E)/(E)dE 

j/(E)dE 


(I) 


where  R(E)  is  the  proton  range  [6].  The  computed  value  is 
<R>  =  4.8  pm  and  most  of  frie  incident  protons  will  stop 
within  4.8  pm  of  surfece,  in  or  near  the  hi^ly  sensitive  gate 
region.  This  distance  corresponds  to  the  range  of  a  300  keV 
proton,  thus  all  protons  with  E  <  300  keV  will  deposit  their 
full  energy  in  the  gate  region,  protons  with  E  =  400  keV  will 
deposit  315  keV  and  the  much  less  numerous  2  MeV 


protons,  130  keV.  Thus,  if  the  external  i^atmis  have  access 
to  the  CCD’s,  the  bulk  of  the  damage  will  occur  in  the  first 
few  microns.  This  is  verified  by  the  work  of  one  groiq>  [7] 
fiiat  has  directly  linked  CTI  increase  to  damage  to  the  buried 
channel  ctmiponent  of  the  CCD,  located  a  fraction  of  a 
micren  beneath  the  gate  region. 

We  have  shown  that  if  the  external  jvotons  have  direct 
access  to  the  CCD’s,  the  peculation  with  energies  in  the 
range  of  100-500  keV  is  reconsiUe  frn*  frie  damage  of  the  FI 
CCD’s.  In  file  next  section  will  sbow  how  the  proton 
transmission  from  enitside  the  spacecraft  to  the  CCD’s  takes 
place. 


in.  Standard  Calculational  Methods 

The  XMM  team  has  used  fiie  Monte  Carlo  code  Geant4  to 
calculate  the  proton  transmissimi  probability.  We  have 
chosen  to  use  to  use  another  well-established  Monte  Carlo 
code,  MCNPX  [8]  to  treat  energetic  proton  scattering.  Both 
Geant4  and  MCNPX  use  cemdensed  collision  physics  to 
compute  the  energy  loss  and  angular  scatteing  of  a  particle 
considering  the  incident  proton’s  cdlisions  with  aUnnic 
electrons  and  with  atomic  nuclei.  Geant4  uses  a  “mixed” 
multiple  scattering  algorithm  [3]  to  predict  {voton  energy 
loss  and  scattering  angle,  x^le  the  MCNPX  physics  for 
determination  of  angular  deflections  is  based  cm  Rossi’s 
Gaussian  model  [9],  and  in  the  enagy  range  of  int^est  here, 
a  continuous-slowing-down  energy  loss  modd.  In  both 
Oeant4  and  MCNPX  models,  more  numerous  collisions  with 
electrons  result  in  small  angular  scattoing  and  a  sonall 
energy  loss  pa  collision.  Less  frequmt  collisions  wifli  tar^ 
nuclei  result  in  comparatively  large  scattering  angles  and 
energy  losses.  Each  collision  in  fids  approadi  is  consido’ed 
to  be  independent  of  the  othras  and  it  is  the  sum  of  all  the 
small,  random  changes  in  angle  and  energy  that  accounts  for 
flie  final  incident  particle  scattering  angle  and  degraded 
energy. 

These  procedures,  and  their  undo'lymg  ]ftiysics  models,  y/ae 
initially  developed  fix  thin  finl  and  bulk  material  scattering 
of  protons  and  alpha  particles  [10].  While  these  scattering 
models  are  valid  fix  moderately  large  angles  of  incidence,  or 
cases  \diere  the  incident  particle  is  traveling  in  the  bulk 
material,  the  fundamental  assumption  of  independence  of  the 
collisions  is  not  valid  fix  grazing  angles  of  incidence.  In 
addition,  the  assumptient  that  flie  incident  particle  will  mta 
into  the  bulk  material,  if  the  particle  trajeeUxy  is  not  exactly 
parallel  with  respect  to  the  surfece,  is  also  not  valid.  In  the 
next  section  we  will  present  a  conceptually  correct  way  to 
handle  grazing  incidence  beam-target  into'actions. 


IV.  Grazing  Angub  Scattering 
A.  Energy  Loss 

In  the  past  decade,  grazing  angle  angular  specular  and  near 
specular  scattering  has  been  studied  as  a  means  of  deducing 
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the  properties  of  the  scattering  sur&ce.  Song  and  Wang  [1 1] 
computed  trajectories  of  grazing  incidence  protons  incident 
(m  a  carbon  sur&ce.  Their  work  ^owed  that,  in  general,  the 
incident  particles  traveled  only  throu^  the  electron  plasma 
cloud  outside  the  sur&ce  before  being  reflected  (no  collisions 
with  atomic  nucld).  Furthermore,  the  lengths  of  trajectories 
inside  the  plasma  cloud  were  only  weakly  dq>endent  on  the 
energy  and  angle  of  incidence.  This  result  provides  an 
explanation  for  the  experimental  results  obtained  by  Winter 
et  al.  [12]  and  Pfonzender  and  Stolzle  [13]  that  the  most 
probable  energy  loss  for  a  wide  range  of  energies  30  <  E  < 
710  keV  is  of  the  order  of  3  keV  and  does  not  dq>end  on  flie 
angle  of  incidence.  In  addition,  the  energy  spectra  of  the 
reflected  protons  are  dominated  1^  a  Gaussian  peak,  centered 
on  the  most  probable  oiergy  loss,  with  only  small 
probabilities  for  greater  energy  losses. 


(aluminum  target).  Data  curve  is  extracted  from  Winter  at 
al.  [12].  Specular  scattering  is  fiir  angle  of  incidence  of  0.5°. 
‘*A11  angles”  is  fiv  scattering  between  0°  and  90°. 

Standard  iQjproadi  to  scattering  calculatimis  cannot 
rq>roduce  the  effect  of  grazing  angle  trajectmies  that  are 
reflected  before  striking  the  surfoce  material.  The 
dif^ences  between  the  two  approaches  are  evident  in  Fig.  3. 
hi  this  figure,  flie  measured  eio’gy  spectrum  of  210  keV 
protons  scattered  off  aluminum  [12]  is  diown  along  with  the 
results  from  an  MCNPX  calculation.  As  expected  the 
MCNPX  results  have  broad,  low  energy  tail  from  the 
numerous  high  energy  loss  collisions  with  the  atomic  nuclei 
in  the  target  material.  A  comparison  of  calculated  and 
measured  energy  losses  for  protons  incident  on  aluminum  as 
a  fimction  of  inddmit  eaeigy  is  shown  in  Fig.  4.  The 
MCNPX  calculated  values  show  a  much  larger  energy  losses 
than  the  data.  The  agreement  is  evm  worse  when  the 
measured  most  probable  values,  which  are  voy  nearly 
average  values,  are  compared  to  MCNPX  average  values. 
The  effect  of  using  the  standard  model  of  scattering  when 
treating  grazing  angle  scattering  is  the  prediction  of  energy 
loss  distributions  fliat  are  highly  asymmetric  and  very  broad. 
This  is  true  for  MCNPX  and  the  fm  the  methods  used  in 
Geant4  [3  ].  The  overall  effect  is  ftw  codes  like  MCNPX  and 
Geant4  to  predict  larger  energy  losses  than  actually  occur  in 
file  scattering  process. 


1000 


S'  100 


<D 


•  Mo«t  Probable 'MCNPX 
o  Avg-ftr  Specular  nafleetion  - MCNPX 
▼  Ay0.1bral«^- MCNPX 
-  Range  of  MaaaurodVWuaa 


1-^ - ^ ^ , - 1 

0  200  400  800  BOO 

Incident  Energy  (koV) 

Fig.  4.  Comparison  of  measured  and  computed  energy  loss  of 
grazing  incidence  protons  (aluminum  target)  as  a  fimction  of 
incident  energy.  Measured  most  probable  energy  loss  values 
[12]  are  all  in  the  region  between  the  two  horizontal  lines 
(all  values  were  qq>roximately  3  keV  in  the  eno'gy  range  of 
50  to  710  keV).  Due  to  the  symmetry  of  the  measured  energy 
loss  distributions  [13],  the  measured  most  probable  value  is 
very  close  to  the  average  value. 


B.  Angulcp’  Scattering 

Angular  scattering  of  the  grazing  incidence  protons  cannot 
be  reproduced  by  physics  models  used  in  Oeant4  and 
MCNPX  These  models  assume  that  the  proton  does  not 
react  wifli  the  scattering  sur&ce  until  it  mters  and  flioi  is 
scattered  as  if  it  were  traversing  bulk  material.  In  &ct,  the 
work  of  Song  and  Wang  [II]  diows  the  incident  proton 
interacts  with  the  sur&ce  long  befiire  it  strikes  it  In  some 
cases  the  proton  is  reflected  befijre  it  strikes  the  surfiice.  For 
larger  angles  of  incidence,  the  proton  only  enters  the  electron 
cloud  that  extends  into  flie  vacuum  to  a  distance  of  a  few 
nanometers  out  from  the  material  surfiice.  In  general,  fiir 
grazing  angle  scattering,  the  incident  protons  do  not  alter 
the  bulk  material  at  all. 

The  angular  scattering  of  grazing  angle  incidence  protons 
from  a  flat  surfiice  is  dominated  by  a  process  vnth  1)  an 
energy  loss  small  compared  to  the  incident  energy  and  2) 
angular  scattering  diaracterized  by  weak  interaction  with  the 
electron  plasma  doud.  In  this  case  flie  proper  scattoing 
response  is  that  given  by  Firsov  [14] 


3 

2tn//  \f/^  +  9^ 


(2) 


where  N  is  the  scattering  response  fimction,  i.e.  particle 
fraction  per  unit  angle  exiting  the  sur&ce,  ^  is  the  angle  of 
incidence  and  0  is  the  scattering  angle.  Note  that  the  results 
are  independent  of  flie  properties  of  bofli  the  scattering 
medium  and  the  incident  particle.  Firsov’s  fiirmula  was 
derived  in  the  limit  of  zero  scattering  energy  loss  and  is 
independent  of  the  precise  fimn  of  scattering  law  as  long  as 
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small  angle  scattering  predominates.  Both  conditions  are 
met  in  the  case  of  grazing  angle  scattering.  It  is  easy  to 
show  that  the  scattering  hinction  N  peaks  at  \{/  =  0  (specular 
reflection)  for  all  angles  of  incidence  and  that  N  is  sharply 
peaked  f<x  small  values  of  \|/  and  q)reads  out  as  \|/  increases. 
Angular  scattering  curves,  calculated  both  using  MCNPX 
and  (2),  for  two  values  of  \\f  are  shown  plotted  in  Fig.  5. 
Given  that  the  vertical  scale  in  Fig.  5  is  logarithmic,  it  is 
evident,  that  the  MCNPX  distributions  are  much  broader  and 
have  a  much  larger  average  scattering  angle  than  the  Firsov 
curves. 


Scattering  Angle  (degrees) 

Fig.  5.  Angular  scattering  distributions  of  120  keV 
protons  incident  on  aluminum  calculated  using  MCNPX  and 
Firsov  distributions  fir  two  angles  of  incidence 

The  consequence  of  MCNPX  predicted  distributions  being 
too  broad  and  being  centered  on  the  wrong  angles  is  that  the 
computation  of  transmission  through  the  telescope  mirrers 
will  lead  to  a  too  low  result  In  reality,  grazing  incidence 
protons  undergo  nearly  specular  reflection  and,  behaving 
very  much  like  x-rays,  are  eflBciently  focused  onto  the  CCD’s 
in  the  fiscal  plane.  The  MCNPX  and  Geant4  scattering 
calculations  predict  that  the  protons  will  be  largely  dispersed 
hy  the  mirrors  and  will  strike  some  other  cmiponent  of  flie 
telescope  and  be  absorbed. 

4.  Conclusions 

This  paper  has  focused  on  scattering  cm  aluminum  rather 
than  the  mirror  materials.  This  is  due  to  the  lack  of 
availability  of  data  for  the  mirror  materials.  However,  the 
overall  grazing  scattering  results  depend  only  weakly  on  the 
scattering  surfiu^e  properties  and  will  not  change 
significantly  for  the  mirror  materials. 

MCNPX  and  Geant4  calculations  underestimate  the  effect  of 
protons  on  focal  plane  instruments  for  two  reasons; 

1)  The  codes  will  overestimate  the  energy  loss  of 
protons,  thus  shifting  to  the  spectrum  of  transmitted 
protons  to  a  higher  incident  energy.  This  effectively 
reduces  the  calculated  number  of  protons  reaching 
the  focal  plane  instruments  since  the 


magnetospheric  incident  protmi  spectrum  in  &lls 
off  rapidly  with  energy, 
and 

2)  The  codes  will  calculate  too  large  scattering  angles 
thus  directly  decreasing  the  number  of  transmitted 
protons. 

The  ideal  solution  would  to  add  the  proper  treatment  of 
grazing  angle  scattering  to  the  MCNPX  and  Oeant4  codes. 
However,  ibis  may  be  a  very  majm-  task  and  not  easily 
accomplished.  A  good  way  of  computing  the  upper  limit  of 
proton  fluence  in  the  focal  plane  is  to  assume  that  the 
protons  undergo  specular  scattering  with  no  energy  loss,  in 
effect  they  behave  like  photons.  This  approach  will  lead  to 
only  a  slight  overestimation  of  focal  plane  flu^ces 
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